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Abstract 

We apply on-shell and integrability methods that have been developed in the context 
of scattering amplitudes in A/" = 4 SYM theory to tree-level form factors of this theory. 
Focussing on the colour-ordered super form factors of the chiral part of the stress-tensor 
multiplet as an example, we show how to systematically construct on-shell diagrams for 
these form factors with the minimal form factor as further building block in addition 
to the three-point amplitudes. Moreover, we obtain analytic representations in terms 
of Grafimannian integrals in spinor helicity, twistor and momentum twistor variables. 
While Yangian invariance is broken by the operator insertion, we find that the form 
factors are eigenstates of the integrable spin-chain transfer matrix built from the mon- 
odromy matrix that yields the Yangian generators. Constructing them via the method 
of R operators allows to introduce deformations that preserve the integrable structure. 
We finally show that the integrable properties extend to minimal tree-level form factors 
of generic composite operators as well as certain leading singularities of their n-point 
loop-level form factors. 
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1 Introduction 

In the last years, there has been tremendous progress in our understanding of AA = 4 Super 
Yang-Mills (SYM) theory in the planar limit primarily based on two different approaches, 
namely the on-shell methods of modern quantum field theory and integrability techniques; 
see [1, 2] and [3] for respective reviews. The former set of ideas and techniques has been suc¬ 
cessfully applied to the perturbative study of on-shell scattering amplitudes of elementary 
states. On the other hand, integrability-based methods, which rely on exploiting all the 
symmetries of the theory, have proven to be very powerful in particular in calculating the 
spectrum of anomalous dimensions of gauge-invariant local composite operators. Though 
some of the integrable structures have also appeared in the study of scattering amplitudes, 
the overlap between both approaches has been rather limited. 

Along with scattering amplitudes and correlation functions, another very interesting 
quantity in a quantum field theory is the form factor, which forms a bridge between the 
previously mentioned on-shell amplitudes and off-shell correlation functions. For a given 
gauge-invariant local composite operator 0{x) in a quantum field theory, the form factor 
is defined as the overlap of the off-shell state created by O from the vacuum |0) at the 
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spacetime point x with an on-shell n-particle state |1,..., n).^ It can be Fourier transformed 
to momentum space, where the operator O carries a momentum q with 7 ^ 0 , yielding 

- ■ ■ ,n]q) = y d^xe"*^'?(l,... ,n|O(x)| 0 ) . ( 1 . 1 ) 

This quantity will be our focus of attention for this paper as it is a perfect candidate to 
study the theory using both the on-shell and integrability techniques. 

In A/" = 4 SYM theory, form factors have been first studied more than thirty years 
back [4], and have received increasing attention of late, both at weak coupling [5-21] and 
at strong coupling [22-24] via the AdS/CFT correspondence. They can be calculated using 
many of the successful on-shell techniques that were developed in the context of amplitudes. 
In particular, BCFW [25, 26] and MHV [27] recursion relations can be applied to construct 
form factors at tree level [5, 7] and the resulting expressions can also be interpreted in 
terms of the volume of polytopes [18]. Form factors have also been studied at loop level 
using generalised unitarity [28-30] not just for the simplest BPS operator and its generali¬ 
sations [5, 11, 17] but also for non-protected operators like the Konishi operator [20], the 
operators in the SU{2) sector [21], and even completely generic operators [19]. All these re¬ 
cent developments have shown that simplicity does exist also for form factors if one studies 
them using the language of modern on-shell techniques. 

However, not all of the interesting features of scattering amplitudes in A^=4 SYM the¬ 
ory have a counterpart for form factors yet. A novel way of studying scattering amplitudes 
has been proposed in [31] using so-called on-shell diagrams, which are bipartite graphs built 
out of two kinds of trivalent on-shell vertices and encode the information of the scattering 
process using fully on-shell data. Moreover, each such scattering process is related to an 
integral over a Grafimannian manifold. In fact, it has been conjectured that all leading 
singularities of the scattering amplitudes of AA = 4 SYM theory as well as the tree-level 
scattering amplitudes can be obtained from an integral on a Grafimannian [32-34]. So far, 
there has been no direct analogue of this geometric picture of scattering for the case of 
form factors. One of the goals of this paper is to provide such a formulation, starting with 
certain tree-level form factors. 

Very little is known about the role of integrability in form factors of A/" = 4 SYM 
theory.^ For amplitudes as well as for correlation functions, integrability manifests itself 
in the appearance of an integrable spin chain at weak coupling. In the spectral problem, 
single-trace operators are mapped to spin-chain eigenstates and the dilatation operator to 
the spin-chain Hamiltonian, see [3] . This integrable Hamiltonian belongs to a whole family 
of commuting operators which also include the corresponding transfer matrices. This family 
can be diagonalised simultaneously using Bethe ansatz methods. More recently, a somewhat 
different spin chain was discovered in the study of amplitudes in A/"=4 SYM theory [35-42]. 
It was pointed out in [43] that the superconformal symmetry and the newly discovered dual 

^As for amplitudes, the on-shell state is specified by the momenta, helicities and flavours of the n 
elementary particles. 

^That is, at least at weak coupling. For an application of integrability to form factors at strong coupling, 
see [23, 24]. 
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superconformal symmetry [44] combine into an infinite-dimensional Yangian symmetry, 
which yields the underlying integrable spin-chain picture [35, 36]. Amplitudes at tree level 
are invariant under this Yangian symmetry. 

Recently, it was proposed that the problem of computing the dilatation operator, i.e. 
the Hamiltonian of the integrable spin chain, in AA = 4 SYM theory can be re-cast in 
a compact form using generalised unitarity methods with form factors being the main 
ingredients [19-21].^ In fact, a special class of form factors called minimal form factors, 
with the number of external fields n equal to the number of fields in the corresponding 
composite operator, realises this spin-chain picture of the spectral problem in the language 
of on-shell super fields used for amplitudes [19]. 

Naively, due to the nature of form factors and also motivated by the results mentioned 
above, we would expect to obtain a relation between the integrable spin chain of the spec¬ 
tral problem and the one that appeared in the study of tree-level scattering amplitudes. 
Indeed, in this paper, we will show that form factors are special states of the latter in¬ 
tegrable spin chain, namely eigenstates of the transfer matrix built from the monodromy 
that yields the Yangian generators studied in the context of amplitudes, provided that the 
corresponding composite operator is an eigenstate of the former integrable spin chain. This 
implies enhanced symmetries for the form factors, analogous to the Yangian symmetry of 
scattering amplitudes. 

This paper is structured as follows. In the remainder of this section, we discuss the 
stress-tensor super multiplet and its super form factors, which we will be studying in most 
of the rest of the paper. In section 2, we briefly review various ideas within the framework 
of on-shell techniques for scattering amplitudes, like on-shell diagrams, their construction 
via BCFW recursion relations and inverse soft limits as well as a Grafimannian integral 
representation, and present the corresponding extensions of these ideas for the case of form 
factors. Our constructions rely on the use of the integrability inspired R-operator techniques 
and the association of a permutation to each on-shell graph as it was done for the scattering 
amplitudes. To allow for a pedagogical presentation, we restrict ourselves to the MHV level 
in this section. Next, in section 3, we further extend the techniques of the previous section 
in order to study similar form factors but at the NMHV and higher N^MHV level. We also 
present some lower-point form factors as examples and conjecture a general Grafimannian 
integral formulation for tree-level form factors. In section 4, we investigate the role of 
integrability for tree-level form factors using the spin-chain monodromy matrix. Specifically, 
we show that all form factors of the chiral stress-tensor multiplet are annihilated by the 
transfer matrix given by its super trace. We also study the action of this transfer matrix 
on the minimal form factors of general operators as well as on on-shell diagrams involving 
them. Finally, in section 5, we conclude with a summary of our results and an outlook 
about future directions. 

Note added On the day of submission, the paper [48] appeared, which has some overlap 
with this article. 

®For the calculation of the dilatation operator from on-shell methods via correlation functions, see 
[13, 20, 45-47]. 
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Form factors of the stress-tensor super multiplet 

In most of this paper, we focus on the form factors of the chiral part of the stress-tensor 
super multiplet, which are the most widely studied ones. Using A/^ = 4 harmonic superspace 
[49] , this part of the stress-tensor super multiplet can be written as 

r(x, 0+) = tr((/)++(/)++) + ■ • • + , (1.2) 

where 9~°‘' = with projectors and . The indices a, a' and 

± correspond to SU{2) x SU{2)' x 1/(1) C 5'[/(4), see [7, 18, 50] for details and further 
references. The lowest component in T{x, 9~^) is the scalar operator with = 

, whereas its highest component is the chiral part of the on-shell Lagrangian 

C. 

The super form factor of this super multiplet is defined as 

• • ,n;q,'y~) = y d'^xd'^(9+“(1 ,... ,n\T{x,9^)\0) , (1.3) 


where 7 “““ is the supermomentum of the multiplet and k denotes the (supersymmetric 
extension of the) MHV degree. For the minimal MHV degree k = 2, the form factor of 
T(x,0+) reads [7]: 


An,2(1, • •. ,n;g, 7 +) 


6\P)6\Q+)6\Q-) 

(12)(23)-- - (n-1 n)(nl) ’ 


(1.4) 


where 

n n n 

P = ^\K-q, Q+ = ^Air)+, Q" = ^ Ai?))" - 7 " (1.5) 

2=1 2=1 2=1 

with ( 5 “'“““ = and 

Note that throughout this paper we will be treating colour-ordered tree-level form 
factors and amplitudes. Hence, we will not indicate this at each expression individually. 


2 The MHV case 

In this section, we demonstrate that many of the recent successful techniques that were de¬ 
veloped for scattering amplitudes can also be applied to MHV form factors, namely on-shell 
diagrams, deformations, R operators and a (deformed) Grafimannian integral representa¬ 
tion. 

2.1 On-shell diagrams, inverse soft limits, BCFW bridges and permutations 
On-shell diagrams 

On-shell diagrams have proven to be a useful tool in the construction of scattering ampli¬ 
tudes. They are built from two different elements, namely the three-point MHV amplitude 

“^The projectors u are related to the m’s by conjugation. 
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^ 3,2 and the three-point MHV amplitude ^ 3 , 1 : 


1 

=^3.2(1, 2 , 3 ) 

1 

=^3.1(1,2, 3 ) 

3 ^ ^2 


-|- A 2 A 2 -l- A3A3)(5®(Ai7yi -|- A2?72 + Asr/s) 

(12)(23)(31) ’ 

( 2 . 1 ) 

^"^(AiAi -|- A 2 A 2 -I- A3A3)i 5^([12] f/3 -|- [23] iyi -|- [31] 7 / 2 ) 

[12] [23] [31] • 


All scattering amplitudes can be built from BCFW recursion relations [25, 26], which can 
be depicted as [31] 



( 2 . 2 ) 


where the BCFW bridge attached at positions 1 and 2 implements the BCFW shift.® Hence, 
they can also be encoded in on-shell diagrams. Similarly to the construction via BCFW 
recursion relations, the on-shell diagram encoding an amplitude is not unique. Equivalent 
on-shell diagrams can be transformed into each other via the so-called square move and 
merge/unmerge move, which are depicted in figure 1 and can be applied to any subdiagram 
of a given on-shell diagram. 









3 4 3 4 



(a) Square move. (b) Merge/unmerge move for black vertices. 

Figure 1: Moves connecting equivalent on-shell diagrams. Similarly to the case for black 
vertices, the merge/unmerge move also exists for white vertices. 

In order to construct form factors via BCFW recursion relations, the minimal form 
factor is required as an additional building block. Hence, it is also required to extend on- 
shell graphs to the construction of form factors. We depict the minimal form factor of T 

®Note that in this work we use the parity flipped version of the BCFW bridge used in [31], i.e. the 
opposite assignment of the black and white vertices in the BCFW bridge. 
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as 


2 


= -^2,2(1, 2 ) 


( 5 ^(AiAi + A2A2 - + A2 ??^)^^(Ai7 /^ +A2^2 ~7 ) 

( 12 )( 21 ) 


(2.3) 

We can then use the construction of the form factors of T via BCFW recursion relations 
[5, 7], which we depict as 



Inverse soft limit 

The MHV form factors J-n ,2 can also be constructed from the minimal form factors via the 
so-called inverse soft limit [51-53] similarly to MHV amplitudes [54]. In total, two types of 
inverse soft limits exist, which either preserve the MHV degree or increase it by one unit. 
In terms of on-shell diagrams, the /c-preserving inverse soft limit amounts to recursively 
adding the structure 

(2.5) 

to two adjacent legs of the diagram. 

For the four-point amplitude ^ 4 , 2 ; this construction starts at the three-point amplitude 
.A 3^2 and can be depicted as 




( 2 . 6 ) 


Similarly, the three-point form factor 
72,2 as® 


1 



73,2 can be constructed from the minimal form factor 



®The result of this construction can easily be seen 
recursion relation (2.4) with a shift in the legs 3 and 1. 


to agree with the 


one obtained from the BCFW 
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Note that the diagram in (2.7) is not invariant under a cyclical relabelling of the on- 
shell legs 1, 2 and 3, although the expression it is encoding is. Hence, we have to add an 
equivalence move for on-shell graphs that involve the minimal form factor, which can also 
be applied to any subgraph of a given on-shell graph. We call it rotation move and depict 
it in figure 2. Together with the equivalence moves for amplitudes shown in figure 1, the 
rotation move guarantees the cyclic invariance of all on-shell diagrams of MHV form factors. 
This is similar to the situation for scattering amplitudes. By itself, also the on-shell diagram 
(2.6) for the four-point MHV amplitude is not cyclically invariant. Its cyclic invariance has 
to be imposed in the form of the square move. However, this suffices to guarantee the cyclic 
invariance of all other MHV amplitudes when combined with the merge/unmerge move. 



2 3 1 


Figure 2: Rotation move for on-shell diagrams involving the minimal form factor. An 
analogous move exists for the inverse combination of black and white vertices. 

Finally, note that we can also construct N™^^MHV amplitudes and form factors via 
the inverse soft limit by adding the ^-increasing structure 

( 2 . 8 ) 

to two adjacent legs of the diagram. The resulting on-shell diagrams are related to those 
of MHV type by exchanging the black and white vertices.' 

Permutations 

For amplitudes, it is possible to associate a permutation 
/I 2 3 ... n \ 

CT= i I 4. ... I = (cr(l),cT(2),cr(3),... ,cj(n)) (2.9) 

V(l) a{2) a{3) ... a{n)/ 

to every on-shell graph by starting at some external particle i and turning right at every 
black vertex and left at every white vertex [31]: 

1 1 

=(2,3,1). (2.10) 

3 " ^'2 3 ^ "' 2 

^In fact, all N''MHV scattering amplitudes and form factors of T can be constructed via the inverse 
soft limit using both (2.5) and (2.8) [54]. In the case of non-extremal k, however, the position and order of 
adding these structures becomes important. 


iT = (3,1,2) 
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If the path ends on particle j, set a{i) = j.® We can extend this to form factor on-shell 
graphs by the prescription to turn back at every minimal form factor: 



( 2 . 11 ) 


As both MHV amplitudes and MHV form factors can be constructed by /c-preserving 
inverse soft limits, we find that their on-shell graphs encode the same permutation, namely 

MHV: cj = (3,...,n,l,2). (2.12) 


Construction via BCFW bridges and permutations 


Using the permutation, a corresponding on-shell graph for general tree-level amplitudes can 
be constructed in a systematic way as follows [31].® First, the permutation is decomposed 
into a chain of transpositions of minimal length, where the multiplication of permutations 
corresponds to the right action, i.e. ai<\a 2 = ((T 2 (cri(l)),..., cr 2 (cri(n))). In this paper, we 
use only adjacent transpositions. Second, each transposition (i, j) is interpreted as a BCFW 
bridge 


j i 


(2.13) 


connecting the legs i and j. Third, these BCFW bridges are applied to an empty diagram 
composed of n lines that start in corresponding vacua,but in the inverse order compared 
to the multiplication in the chain of transpositions; this means that the rightmost trans¬ 
position corresponds to the BCFW bridge that is applied first to the vacua.Fourth, the 
vacua, the edges starting at the vacua and every vertex that is connected to less then three 
edges is removed to obtain an on-shell diagram. This construction is illustrated for A 3,2 in 
figure 3. 



Figure 3: Constructing the on-shell diagram of A 3,2 via permutations and BCFW bridges. 

®In contrast to [31], we are using ordinary permutations instead of decorated permutations here; see the 
discussion below. 

®We use the graphical notation of [40]. 

^°We will give further meaning to these vacua below; they are the same as the “lollipop” diagrams of [31]. 
Note that the kind of the different vacua, i.e. -I- or —, is imposed by hand following the prescription of [40] 
and not given by the permutations, which are not decorated in this work. 

^^The inverse order follows from our choice of BCFW bridge. 
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The on-shell diagrams of the MHV form factors can be systematically constructed via 
BCFW bridges using essentially the same construction as in the amplitude case. The only 
difference to the amplitude case, where only one-site amplitude vacua appear, is that the 
minimal form factor occurs as a vacuum at positions n and n — 1. This construction is 
illustrated in figures 4, 5 and 6. 


(3,1,2) = (2,3) <(1,2) 

3 2 1 2 

Figure 4: Constructing the on-shell diagram of J -^^2 via permutations and BCFW bridges. 




Figure 5: Constructing the on-shell diagram of J- 4^^2 via permutations and BCFW bridges. 



Figure 6: Constructing the on-shell diagram of J- 5^2 via permutations and BCFW bridges. 


2.2 Deformed form factors and R operators 

We can now introduce deformations of the form factors and construct these deformed form 
factors in analogy to the amplitude case [35, 37, 39-42, 55]. For amplitudes, a sequence 
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of BCFW bridges can be translated into a chain of R operators that acting on a suitable 
vacuum state produce a deformed version of the amplitude, or rather some BCFW term 
of it. In this section, we will use the R operators primarily as means to obtain analytic 
expressions for the form factors, in particular representations in terms of Grafimannian 
integrals. However, the R-operator formalism is based on the spin-chain picture of integra- 
bility and we will use this fact in section 4 to show that form factors are well defined states 
in the integrable model and posses enhanced symmetry properties. There, we will also give 
further details concerning the definition and the properties of the R operators and discuss 
the integrability-preserving deformations. 

The R operators [35] can be defined by their action on general functions / of the 
kinematic data,^^ 

/ ci C-^ 

- a\j,Xi,fn,Xj,Xj + aXi,fij + arij) . (2.14) 

Here, the parameter u will eventually correspond to a (integrability-preserving) deformation 
of the physical form factor. Moreover, the vacua that occurred in the previous discussion 
are given by 

+ ~ 

= 5+ = <I2(A0, I = 5-= 6\Xi)6\vi). (2.15) 

i i 


There are two types of them, reflecting the different possible MHV degrees of the final 
expression. 

Let us consider the three-particle MHV amplitude >13^2 in figure 3 as example. The 
sequence of transpositions (2,3) 0 (1,2) translates into 


R23('W32)Rl2(lt3l)<^i^<^2 <^3 


(12)1-“23 (23)1-“3i (31)1-“i2 


(2.16) 


where Ui are parameters associated to deformations of the local central charges (see section 
4) and 

Uij = Ui — Uj . (2-17) 

The undeformed three-particle MHV amplitude A 3^2 is recovered in the limit Ui —>■ 0. 

The previous discussion suggests that we can use essentially the same construction for 
the three-particle MHV form factor; the only diagrammatic difference is the substitution 
of the minimal form factor for the vacua at sites 2 and 3, cf. figure 4. Using the minimal 
form factor (2.3), which we label by the two sites it occupies, we find that the same chain 
of R operators produces a deformed version of the three-point MHV form factor 


R23('W32)Rl2(w3l)<^l’-72,2(2, 3) = 


- q)6\j:U - 7“) 


(12)1-“23 (23)1-'“3i (31)1-“i2 


(2.18) 


In the limit of vanishing deformation parameters, this reduces to (1.4) with n = 3. 


^^Note that we extend the usual definition to harmonic superspace. 

^®Here and in what follows, we will ignore phases that also appear in the amplitude case. 
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Since all n-point MHV form factors can be obtained by iterated inverse soft limits, this 
construction generalises to all n, in particular to the further examples shown in hgures 5 
and 6. The result is 


... ,n) = 


AA, - A.r?+), 5 "(Er=i A»r?- - 7") 

nr=i(* i+i)i-“*+ii+2 


(2.19) 


which reduces to (1.4) in the limit of vanishing deformation parameters. 

Instead of performing the construction of (2.19) via R operators explicitly, we will now 
use the R operators to obtain a Grafimannian integral representation that is valid for all 
n and evaluates to (2.19). 


2.3 A (deformed) Grafimannian integral representation for the MHV form 
factor 

The minimal form factor (2.3) can be rewritten in a form that closely resembles the vacua 
(2.15) used in the construction via R operators: 

7 - 2 ,2(1,2) = = 6(2 ■ (2.20) 


Here, we have absorbed the off-shell (super) momentum of the operator into modified 
kinematic variables for the on-shell states. 


(2k 




(21)’ 

2i - hi ^21) ’ 

2i 

- 2i 

(Ik 

-- -- (1|7" 



(12)’ 

22 -22 ^12) ’ 

22 

- 22 


( 2 . 21 ) 


Note that this expression looks exactly like 61 5^, though with twisted kinematics A and ^ 
that contain the information about the off-shell (super) momentum insertion. 

Using this form of the minimal form factor, we can apply the same sequence of R 
operators as in (2.18) to obtain, before integration. 


•^3,2(1 , 2, 3) — R23('*^32)Ri2(«3i)<^i '^23 

= J J+U32 J ^^iC{ai,a2) - A) ^^(C'(ai,a2)-2) 5^{C-^{ai,a2) ■ X), 

^ ^ ( 2 . 22 ) 
where the super spinor helicity variables 2 and 3 are twisted as in (2.20) while 1 is untwisted. 
The matrices C and C"*" are orthogonal to each other, i.e. C{C'^)'^ = 0, and given by 


C'(ai,a2) “ ^ Q, ’ C^{oiit012) — (l —ai 0102 ) ■ 

Their products with the external super spinor helicity variables are defined as 

{c ■ A)? = E Cult , {c-n)f = jz Cunt , • A)} = E cjx , 

2=1 2=1 2=1 


(2.23) 


(2.24) 
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where / = 1,..., fe and J = 1,..., n — fc. We can also write (2.22) in a GL{2) invariant way, 
as an integral over the Grafimannian G(2,3): 


-^3,2(1, 2, 3) 


d2x3(^ 

(12)1-“23(23)1-«3i(31)1-«i2 




(2.25) 


where (ij) denotes the minor of C that is built from the columns i and j. This is precisely 
the (deformed) Grafimannian integral for the three-point MHV amplitude [41, 42] with the 
twisted kinematics accounting for the operator insertion. 

We can generalise the above derivation to an arbitrary number of external on-shell 
fields: 


7'„,2(1, ... ,n) 


nr=i(* i + 


6 ^{C-X) 6^{C-u) 


(2.26) 


Here, the shifted kinematic variables can actually be at any two positions. One can easily 
check that the sequence of R operators necessary to derive this expression does not contain 
BGFW shifts that would spoil this simple dependence on the modified kinematics A, fi. It is 
also trivial to check that this integral representation gives the correct result upon localising 
the integration on the support of the delta functions: we simply take the (deformed) Parke- 
Taylor formula and replace the kinematic variables, A ^ A, r) —)• 2. Since the A’s are 
not modified, the Park-Taylor prefactor is unaffected by this replacement, and the only 
effects are shifts in the (super) momentum conserving delta functions, P ^ P — q and 
Q~ —>■ Q~ — 7“. This follows from the identity 




+ A“ 



\Ci\Q: I \Ci\a 


e,0{XfX]-X-X]) 

{ji) 


/3 a 


=5% 


(2.27) 


and a similar identity for the 77’s. The above argument shows that (2.26) correctly repro¬ 
duces (2.19) and in particular the undeformed result (1.4). 


3 Beyond MHV 

In the previous section, we have considered the simplest form factors, namely the MHV 
form factors, to introduce many important concepts. In this section, we will see that these 
concepts continue to apply beyond MHV, although with some modifications. In particular, 
we will conjecture a Grafimannian integral representation for all form factors, both in spinor 
helicity as well as in twistor and momentum twistor form, and provide several non-trivial 
checks. 

3.1 On-shell diagrams and R operators 

Since all form factors can be constructed via BGFW recursion relations as shown in (2.4), we 
can also directly associate on-shell diagrams to each BGFW term — completely independent 
of the MHV degree k. One main difference between k = 2 and k > 2 is that all MHV 
form factors can be constructed via the inverse soft limit without regard to the order and 


13 








insertion positions, which directly gives the on-shell diagram. For k > 2, the result of the 
BCFW construction, and hence the correct on-shell diagram, is less obvious. A second main 
difference between k = 2 and A; > 2 is that, both for amplitudes as well as for form factors, 
there are several BCFW terms and hence on-shell diagrams which have to be summed to 
obtain the complete expression. However, for amplitudes, they can be combined into a single 
top-cell diagram, which corresponds to a top-dimensional integral over the Grafimannian 
and yields all required BCFW terms when taking suitable residues. We will find in section 
3.2 that we can define such top-dimensional integrals also for form factors. However, we 
will see below that a sum of several top-cell diagrams will be required. 

Likewise, it is always possible to construct a given on-shell graph by acting with a chain 
of BCFW bridges on suitable vacua. Translating these BCFW bridges to R operators, we 
can build deformed BCFW terms and top-cell diagrams for form factors as in section 2.2. 
Hence, in order to construct Grafimannian integrals and R operator representations, the 
first important step is to identify the corresponding top-cell diagrams. Let us look at several 
examples first. 

jVjmaxMHV 

A special class of form factors beyond k = 2 \s given by N™*^^MHV, which has k = n. 
For amplitudes, the corresponding case is MHV, which has k = n — 2. Similarly to MHV 
amplitudes, N™'^^MHV form factors can be constructed via the inverse soft limit without 
regard to the order and insertion points. Hence, the on-shell graph is immediate. The 
permutation associated to these on-shell diagrams is 


NmaxMHV : a = (n-l,re,l,2,...,n-2). 


(3.1) 


In the construction via R operators, we now have n — 2 conjugate amplitude vacua 6 ^ on 
the right of the minimal form factor instead of n — 2 amplitude vacua 5f on its left. In the 
simplest case of n = A: = 3, this is depicted in figure 7. 



Figure 7: Constructing the on-shell diagram of via permutations and BCFW bridges. 


NMHV 

The first case that is truly beyond MHV is A: = 3 for re > 4. For the case of re = 4, 
the BCFW sums for all adjacent shifts are shown in figure 8. They have been generated 
using (2.4). Applying the moves in figures 1 and 2, it is easy to see that 

^{i+2) mod 4 ) ^(i+2) mod 4 • (3-^) 
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Figure 8: BCFW terms of ^”4^3 for all adjacent shift. The line stems from a shift in i 
and i + 1 . 


These BCFW terms can be obtained as residues from the sum of two different top-cell 
diagrams. The first of these is shown in figure 9 together with its permutation and con¬ 
struction via R operators; the second one can be obtained from it by a cyclic shift of the 
external on-shell legs by two. Concretely, all vertical edges in the top-cell diagram in figure 
9 are removable. Deleting them, we obtain from left to right C2, Ai, C 3 and ^2- 

Several remarks are in order. First, we do require more than one top-cell diagram to 
generate all BCFW terms. Second, the top-cell diagram is not cyclically invariant, and nei¬ 
ther is the corresponding permutation. Instead, we (in principle) have to explicitly consider 
all cyclic permutations of the top-cell diagram and the corresponding permutation. Third, 
the permutation is not decomposed into a minimal number of transpositions. 

It is possible to construct other cases with higher n, k in an analogous way. 

N^MHV and a relation to amplitude on-shell diagrams 

We conclude this subsection with a general observation relating the on-shell diagrams of 
form factors with those of amplitudes. In particular, this will lead to (a conjecture for) the 
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(4, 2,3,1) = (1,2) < (3,4) < (2,3) < (1,2) < (3,4) 



Figure 9: Constructing the top-cell diagram of via permutations and BCFW bridges. 
Note that the decomposition of the permutation that produces the correct on-shell diagram 
is not minimal in the amplitude sense. 

form factor top-cell diagrams at general n, k. 

We note that the n-point form factor shares interesting features with the (n-|-2)-point 
amplitude. To begin with, the MHV degree k ranges from 2 to n in both cases. Moreover, 
n -|- 2 is the expected number of kinematic dependencies if we consider that the off-shell 
(super) momentum of the operator can be parametrised by two on-shell (super) momenta. 

Finally, we have found that we can obtain the top-cell diagrams of J-n,k from the top¬ 
cell diagram of An+ 2 ,k by applying moves until a box appears and replacing this box with 
the minimal form factor. Graphically, this relation reads 

n -I- 1 n -I- 2 


(3.3) 


^•••321 n---321 

where we have replaced the box at the legs re -|- 1 and re -|- 2 for the sake of concreteness. 
This relation is valid for all form factors presented in this paper. 

At the level of the BCFW terms, a respective relation can be seen to be true by con¬ 
sidering the BCFW recursion relations (2.2) and (2.4). Note that for re = 2 the amplitude 
on-shell diagram is nothing but a box and entirely replaced by the minimal form factor. 
Constructing the amplitude An+ 2 ,k recursively via (2.2), we find that boxes can only occur 
at the boundary of the on-shell diagram. Then constructing J-n,k recursively via (2.4), we 
find that each term in (2.4) can be obtained by replacing one of the boxes in a term of the 
construction of An+ 2 ,k by the minimal form factor. It would be interesting to prove this 
relation also at the level of the top-cell diagram(s). 

Using the relation (3.3) between Tn,k and An^ 2 ,k, we can also relate the corresponding 
permutations. By replacing the box by the minimal form factor, we can hide the corre¬ 
sponding legs of the amplitude in the composite operator. At the level of permutations, 
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this connects the preimage of the hidden leg to its image. For ^n+2,A:i the permutation 
corresponding to the top-cell diagram reads 

A+2,fc : a = {k + l,... ,n,n + l,n + 2,l,2,... ,k) . (3.4) 

If we hide the legs n -|- 1 and n +2, we obtain 

^n,k'- cr = (A:-I-1,... ,n, A: - 1, fe, 1,2,..., - 2). (3.5) 

Moreover, the top-cell diagram with the minimal form factor replaced by the open legs 
n -|- 1 and n -|- 2 is characterised by 

J^n,k without J'2,2 : d = (A: -|- 1,..., n, n -|- 2, n -|- 1,1,2,... , A:, A; — 1). (3-6) 

The permutation a allows us to directly generate this on-shell diagram e.g. using the 
Mathematica package positroid.m [56]. 

Let us conclude with a comment about the role of the permutation for form factors. 
As for scattering amplitudes, it is invariant under all equivalence moves. In contrast to 
the case for tree-level scattering amplitudes, the permutation a for the complete form 
factor on-shell diagram requires a decomposition into more than the minimal number of 
transpositions to construct the corresponding on-shell diagram via BCFW bridges as shown 
above.However, the modified permutation a can be directly used to obtain the on-shell 
diagram with the minimal form factor replaced by the open legs n -|- 1 and n -|- 2. This is 
similar to the situation for one-loop amplitudes, whose on-shell diagrams are also best not 
constructed from their permutations but from the permutations of the corresponding higher 
point amplitudes before taking the forward limit. It would be interesting to fully classify 
the equivalence classes of form factor on-shell diagrams and find appropriate combinatorial 
labels for them. 

3.2 Grafimannian integrals for higher MHV degree 

Having identified (a conjecture for) the general top-cell diagram, let us now turn to the 
Grafimannian integral representation for form factors. 

General considerations on the Grafimannian 

The fundamental idea behind Grafimannian integral representations of scattering ampli¬ 
tudes [31-33] is to express momentum conservation in a geometric way.^^ Regarding the 
external kinematic data as a pair of two-planes A and A in n-dimensional space, momentum 
conservation is expressed as the orthogonality of these planes: 

n 

= = (3.7) 

i=l 

would be interesting to see whether the decomposition is minimal when adding further conditions 
such as considering 1 and n to be non-adjacent for the purpose of the decomposition. 

^®Here, we focus on momentum for brevity. The same arguments apply to super momentum by replacing 
A with fj. 
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The Grafimannian representation linearises this constraint by introducing an auxiliary hy¬ 
perplane C G G{k, n) such that 

n n 

{C-X)j = Y,CuK = ^ and • A)} = = 0 ^ A-A = 0, (3.8) 

i=l i=l 


where C*-*- is the orthogonal complement of C fulfilling C'(C'“*“)^ = 0 and / = 1,..., A:, J = 
1,..., n — A:. The Grafimannian integral for scattering amplitudes integrates a holomorphic 
form on G{k, n) on the support of these constraints. 

As discussed in section 2.3, we can similarly geometrise momentum conservation for 
form factors. Setting 


Afc — ) k — 1, 


, n. 


k^i,j 


Az — 


01 g 

00 


Aj — 


00 

( 0 ) 


7^ (3-9) 


for arbitrary i and j, we can express momentum conservation as A • A = 0; cf. (2.27). 
We also saw that in the MHV case a naive way of introducing an auxiliary Grafimannian 
works. One can simply use the same Grafimannian and the same form as one would use for 
the MHV amplitude with the same number of legs. It is clear, however, that this way of 
linearising the geometrical constraint cannot work beyond MHV. For instance, the MHV 
degree k ranges up to n for form factors, while for amplitudes it only ranges up to n — 2. 
Since G{n, n) is just a point, a larger Grafimannian is necessary for N™®'^MHV, but in fact 
already starting from NMHV. 

The different range of MHV degrees already suggests that the correct Grafimannian is 
G{k,n + 2); this also fits nicely with the observation (3.3) as well as with the general fact 
that an off-shell momentum can be parametrised by two on-shell ones. Indeed, instead of 
(3.9), we can define new kinematic variables as a pair of two-planes in an (n-|-2)-dimensional 
space as 

Ai — ) 

Ai — ) 

where and are arbitrary non-collinear reference spinors. Momentum conservation is 
then expressed as A" A = 0. As shown in (2.27), the two additional on-shell momenta indeed 
encode the off-shell momentum: An-i-iAn-i-i + An-i- 2 An-i -2 = —g- Analogously, we can define 
fermionic variables 


f = 1,.. 

■ ,n , An-l-l — ^A 1 

An+2 

f = 1,.. 

\ (^s|g 

An+2 


{U\q. 
Ha^b )’ 


(3.10) 


ut = vt > 
m = > 


i = 1,.. 
i = l,.. 


n , ut+i = 0 > 


n, Un+l = 


{^bCa) ’ 


Sn+2 = 0 > 

-- _ _(Mzl 

(U^b) 


(3.11) 


which encode the off-shell super momentum as Xn+iILn+i + hn+2Un+2 — ~1~ ■ Super¬ 
momentum conservation can then be written as A • 2 = 0. 

We can now linearise the constraint imposed by (super) momentum conservation by 
requiring C" • A = 0, C" • g = 0 and G'^ • A = 0 with G' G G{k, n + 2). 
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Prom on-shell graphs to Grafimannian integrals 

To show that form factors can be written as integrals over the Grafimannian G{k,n + 2), 
we break the corresponding diagram into two pieces: the minimal form factor (2.20) and 
a purely on-shell piece with re -|- 2 legs for which a Grafimannian integral representation is 
known. We then glue these two pieces together, i.e. we perform the on-shell phase space 
integration. We start by discussing this procedure in a general form that can be applied to 
any diagram; then, we will look at some low-point examples to see how the explicit form 
for top-cell diagrams looks like. 

The on-shell piece that will be glued with the minimal form factor can be written 
as [31] 

/ = / ^ ... ^ 5^x2((7 . A) ■ fj) 5(^+2-fc)x2^(^± . ^ 

J (X\ 

where the matrix C depends on the ccj’s, C = C{ai) € G{k, n + 2) and m is the dimension 
of the corresponding cell in the Grafimannian. Gluing the minimal form factor to the legs 
re -|- 1 and re -|- 2 corresponds to calculating 



(vol(OL(l)l ’'7 




T 

n+l n+2 


A^-A 


1(1,..., re -|- 2), 


(3.13) 


where was defined in (2.20) and the signs of the corresponding A’s are inverted since 
the two particles are ingoing with respect to . We can perform the A and i) integration 
using the delta functions of the minimal form factor (2.20); this replaces 




^n+2 


(re -7 2 \q 
{n+2 n+l) ’ 
(re -I- l\q 
(re-|-l re-|-2) ’ 


Vn+l 


Vn+2 


(re -h 2|7 
(re-|-2 re-|-l) ’ 
(re -h 1|7“ 
(re-|-l re-|-2) ’ 


vt+i 0 , 

Vn+2 0 • 


(3.14) 


To remove the GL{1)‘^ redundancy in the remaining A integrations, we parametrise 


An-l-l = ^A- (blis , A„+2 = is - hiA , 


(3.15) 


where and + are two arbitrary but linearly independent reference spinors, which will 
be identified with the ones in (3.10). With this, (re-|-l re-|-2) = (/li/32 — l)(?B'^2l) and the 
replacement (3.14) becomes 


Vn+l 

^n+2 

Vn+2 


1 


_ {CbIq 

A/32 - 1 {(b(a) 
1 {^Bh~ 

A/32 - 1 {(bCa) 

1 {U\q 

A/32 - 1 (?a?b) 
1 (?a|7’ 


P2 {U\q 


/3+2 - 1 (Ca^b) 

h (^a|7~ 
Pih-l {^a(b) 
Pi {^B\q 
P1P2 - 1 {^bCa) 

Pi (^b|7~ 

P1P2 - 1 {^aCb) ' P1P2 - 1 {^b^a) 


(3.16) 


+ 
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At the same time, the measure transforms to 


d^A„+i 


d^A- 


nH-2 


= {CaCb){Cb^a) / d^id/32 


J Vol[GL(l)] Vol[GL(l)] 

Applying the substitutions in (3.16), we can write (3.13) as 

dam d/ll d/32 


(3.17) 


It = {^a^b){^bCa) f —- 

J ai 


Om (1 - /3i/32)2 

X d^^\C'{ai,/3i) • A) 6’^^\C'{ai,/3i) • § /3i) • A) , (3.18) 

where we recombined the columns of the matrix C snch that they form the coefficients of 
the kinematic data A, U defined in (3.10) and (3.11). This new matrix C = (CJ • • • C'^_^ 2 ) 
depends both on the aj’s as well as the /3j’s. Its hrst n columns coincide with those of C 
and the last two columns are given respectively by 


1-/3i^2 


Cn+l + 




1 




1 — /dl/32 


<^^+2 , C^+2 ^ Cn+2 + ^ Cn+l ■ 

(3.19) 

Hence, also the first n columns of C'-^ coincide with those of C-*- and the last two columns 
are 

C;ii = Ci+i-/32Ci+2, =C^+2-AC^n+i- (3-20) 

The factor of (1 —/3i/32)^ in (3.18) is a Jacobian from reorganising the C-^-X delta functions, 
which we can write as 

^(n+2-fc)x2(^± • A) = n / A(’^+2 )x2 ^ (3 21) 

K=l'^ 

where with AT = 1,..., fe is a set of auxiliary variables. For the delta functions corre¬ 
sponding to the columns n + 1 and n + 2, we have 

<3^(An-|-l — PlCl n+l)S‘^{^n+2 — PlCl n+ 2 ) 


A^(An-Hl - PlC'l n+1 - /3l(An-H2 “ PlC'l n+ 2 )) 
^^(An-l-2 — PlC'i n+2 ~ ^2(An-|-l “ PlC'^ „+i)) 
1 


(3.22) 


(1 - /3i/32)2 


5 (An-l-1 — PlC'i „+i) j (An-l-2 — PlC'l n+ 2 ) 


The above shows that diagrams contributing to the form factor can be expressed natu¬ 
rally as some Grafimannian integrals. Of course, it remains to find some general expression 
for the form that is to be integrated over the support of the delta functions. For this we will 
look at some concrete examples first; although the gluing procedure works for any on-shell 
diagram, we will from now on focus on top-cell diagrams based on the conjectured relation 
with the amplitude diagrams in section 3.1. 
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MHV revisited 


In this section, we will show that the general procedure outlined in the last section reduces 
to the results from section 2 for MHV degree k = 2. This will also give us a first idea of 
how the form of the G{k, n + 2) Grafimannian integral looks like in the general case. 
Consider the C' matrix in the standard gauge fixing, 


A 0 c ']^3 • • • 

\0 1 C23 ■ • • C2 n +2 j 


(3.23) 


and kinematic data A, A, ^ with the off-shell information encoded at positions n + 1 and 
n+2. Compared to the G{2, n) integral given in (2.26), we have four additional integrations 
as well as four additional bosonic delta functions 5{G''^ ■ A) involving the (n + 1)**^ and 
(n + 2)*'^ rows of C. If we choose the reference spinors is such that = An-i-i = A 2 
and = Xn +2 = , these four additional delta functions impose 


- cA+i A? - Af + A^ = 0, -cA+2 a? - c'^n+2 Af + A? = 0 . 


(3.24) 


Upon integrating out these delta functions, we obtain a Jacobian (12) ^ which cancels the 
prefactor in the general expression (3.18), and the G' matrix is set to 


AOc '3 c'„0l\ 

i^oic' 3 --- cAioy 


(3.25) 


By defining G as G' without the last two columns, the delta functions are now identical to 
the ones in (2.26), 

A(C-A)5®(C-A ■ A) (3.26) 

with the twisted kinematics at position 1 and 2, as in (2.21). 

We now have to check whether the form obtained by the gluing procedure yields the 
same form as in section 2 after this integration. We have performed the gluing explicitly for 


I = 



d = (3,5,4,2,l). 


C 


C l 0 —Q ;2 —0(204 —Oi 

0 1 03 0304 0 


Figure 10: On-shell sub-diagram I obtained by removing the minimal form factor from 
the on-shell diagram of shown in figure 4, corresponding permutation a and G matrix. 

the MHV form factors with up to six external particles. The on-shell sub diagrams obtained 
by removing the minimal form factor from the on-shell diagrams, the corresponding permu¬ 
tations a as well as the G matrices obtained from the Mathematica package positroid.m 
[56] are shown in figures 10 and 11 for 3 and 4 points, respectively. We invariably found 
that after changing from edge variables a,, /3i, (32 to canonically gauge-fixed Cij variables, 
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0 —03 —02 — O3O5 -05(02 + O3O5) —Oi 
10 1 04 04O5 04O5O6 0 


Figure 11: On-shell sub-diagram I obtained by removing the minimal form factor from 
the on-shell diagram of J^4^2 shown in figure 5, corresponding permutation a and C matrix. 


the integral could be written in the following form:^® 


where 


^2in+2)fjf Y{1 - y)-l 

Vol[GL(2)] (12)(23)---(n+ln+2)(n+2 1) 


(3.27) 


y _{n n+l){n+2 1) 

(n n-|-2)(n-|-l 1) ’ 


(3.28) 


We checked up to n = 6 that, after integrating out the four additional delta functions as 
in (3.24), the form reduces to one from section 2.3: 


y(i -y)-i 

(12)(23) ■■■{n n+l){n+l n+2){n+2 1) 


1 

(12)(23)-- - (n-1 n)(nl) 


(3.29) 


Three-point NMHV 


I = 



d = (5,4,l,3,2). 


C = 


^ 1 03 02 0 0\ 
0 0 1 oi 0 
ya4 0 0 0 1/ 


Figure 12: On-shell sub-diagram I obtained by removing the minimal form factor from 
the on-shell diagram of A'3,3 shown in figure 7, corresponding permutation a and C matrix. 

The simplest NMHV form factor is ^3,3. Diagrammatically, it can be obtained from a k- 
increasing inverse soft limit of the minimal form factor. Using the general gluing procedure 

^®Note that we have ignored overall signs in the gluing procedure since the sign of residues expressed in 
edge variables is not readily determined. See [57] for an elaborate algorithm that determines these signs, 
though. 
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outlined above, we find after a change of variables the following GraBmannian integral 
representation: 




d 3 x 5 (^/ y(i-y)-i 

Vol[GL(3)] (123)(234)(345)(451)(512) 


6 \C' -X) 


(3.30) 


where 


(234) (512) 

(235) (412) 


(3.31) 


and the off-shell (super) momentum is encoded in the on-shell variables at position 4 and 
5; cf. (3.10). 

Let us now evaluate the Grafimann integral (3.30). After gauge fixing, the matrix C' 
reads 

/I 0 0 c'lgX 

C' = 0 1 0 c'24 ^5 • (3-32) 

Vo 0 1 c'34 cU 


We can solve for 0(4 and c '5 with i = 1,2,3 by contracting the terms inside 6 ^{C' ■ A) with 
A 4 and As. This yields 


r' = -M = i^B\q\i] 

[ 45 ] g 2 ’ G5 ^2 


(3.33) 


where we have used (3.10) in the second step. Inserting (3.33) into 5^(G'-*--A), we obtain the 
momentum-conserving delta function contracted with A 4 and A 5 . Undoing this contraction 
yields a Jacobian of [45]^, which, together with the Jacobian [45]from the previous 
contraction with A 4 and A 5 , gives [45]”^.^^ 

Inserting the solutions (3.33) into (3.30) and applying the Schouten identity, we find 


{CaCb)‘^ 

(q 


d3x5(^/ 


y(i-y) 


-1 


Vol[GL(3)] (123)(234)(345)(451)(512) 

J^2(C"-2) S'^C^Pi-q), 


d^(C'-A) d^^(C'-u) d^(C'^-A) 


(3.34) 


[12] [23] [31] 


i=l 


which agrees with the result of [7]. 

Note that the cyclic invariance of the form factor is not manifest in (3.30). The final 
expression (3.34) obtained from its evaluation, however, is manifestly invariant under cyclic 
relabelling of the legs 1, 2 and 3, as can be seen from (3.33). 


Four-point NMHV 

As discussed in subsection 3.1, the four-point NMHV form factor is the first example of a 
form factor for which it appears natural to combine different BGFW terms diagrammati- 
cally into a top-cell diagram with additional edges. Note that the gluing procedure outlined 
in the beginning of this subsection together with the connection (3.3) will indeed generally 
lead to a top-dimensional integral over the GraBmannian G{k,n + 2). 

^^Note that we have dropped the double underscore in the notation for the spinor brackets. 
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4 ^ 1 




I = 


-O- 


a = (4, 6 ,5,1,3,2), 




C = 


(1 0.2 ~\~ CH4 0 — O 2 O 3 — 020 ^ 0 ^ 0 

0 10 —03 —o^oq —oi I . 

\0 0 1 05 + oj o^oq 0 


Figure 13: On-shell sub-diagram I obtained by removing the minimal form factor from 
the top-cell diagram of ^^ 4^3 shown in figure 9, corresponding permutation a and C matrix. 

From the general expression (3.18) and the on-shell diagram and C matrix shown in 
figure 13, we obtain a result that can be written in the following form: 

. ^3x6(27/ 


where 


{(aCb)" 


^ 4,3 — 


Vol[GL(3)] 


04,3 5 ^{ C ' ■ A) • 2) ■ A), 

(345)(612) 


Y = 


(3.35) 


(3.36) 


(123)(234)(345)(456)(561)(612) (346)(512) 

Gluing the same diagram at legs 2 and 3 and relabelling to obtain the other top-cell 
diagram, we find with the same A) A) U for the first diagram: 




Vol[GL(3)] 




4,3 


1 2 3 4 5 

4' 'I' 4' 'i' 4' 'I' 

3 4 12 5 


5®(C' • A) ■ u) 5“(C'^ • A) • (3.37) 


srl 2 /^/ 




Conjecture for all n, k 

Based on the (conjectured) relation (3.3) between the top-cell diagrams of amplitudes and 
form factors for generic re, we have computed the Grafimannian integrals for all form 
factor top-cell diagrams up to six points, with the minimal form factor glued at positions 
re -|- 1 and re -|- 2. The on-shell diagrams that need to be glued in this case are labelled by 
the permutation given in (3.6). We have invariably found the following representation: 

P ^kx{n+2)(jl 


{^AisY 


\ol[GL{k)] 


^n,k ■ A) ■ u) 52 x(n+ 2 -fc)(^/± . ^ 


(3.38) 


where 


^n,k — 


Y = 


Y{l-Y) 


-1 


{1 ■ ■ ■ k){2 ■ ■ ■ k+1) ■ ■ ■ " k — 3){n+l ■ ■ ■ k — 2){n+2 ■ ■ ■ k — 1) ’ 


{n — k+2 ---n n+l){n+2 I - ■ ■ k — l) 
{n — k+2 ---n n+2){n+l I - ■ ■ k — l) 


(3.39) 


Above, we have explicitly checked that this conjecture leads to the correct result for k — 2, k = n = 3 
and k = n — 1 = 3. Further checks will be given below. 
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and the off-shell data is encoded in the kinematical variables at the position n-|-1 and n + 2 
as in (3.10).^® In general, up to n copies of this form have to be considered, which arise 
from cyclic shifts in the labels 1 to n. 

We will present further checks for this conjecture in the next subsection. 

Note on deformations In contrast to the MHV form factors considered in section 2, we 
did not employ the method of R operators for more general form factors and hence we have 
so far not included deformations in the Grafimannian for N^MHV as counterpart to what 
was done for amplitudes in [41, 42]. The reason we preferred a direct gluing procedure is 
that it immediately leads to a top-dimensional integral over the Grafimannian. The method 
of R operators, while still applicable, will in general result in a Grafimannian integral with 
some of the delta functions already integrated out. Nevertheless, a general expression for 
the deformed top-dimensional form could be obtained in this way; we leave this for future 
work. 


3.3 Twistor and momentum twistor Grafimannians 

Next, let us transform the previously obtained Grafimannian integral representation from 
spinor helicity to twistor and momentum twistor variables. Throughout this subsection, all 
kinematic variables are defined as in (3.10) and (3.11). In order to facilitate notation, we 
will hence drop the double underscore from all spinor (and twistor) brackets. 


Twistor space 

Given the Grafimannian integral (3.38) in momentum space, we can transform it to twistor 
space in analogy to what was done in the amplitude case in [32]. 

The super twistor space we use here corresponds to our special choice of spinor helicity 
variables (3.10) and (3.11) and is given by the set of all super twistors 

yVi = , (3-40) 


where is related to A* via Witten’s half Fourier transformation [58] 


J exp(—• . 


Via (3.41), the prefactor in (3.38) can be written as 


= { 


d 


d 


d^n+l d^n+2 


r- 


(3.41) 


(3.42) 


The delta function 52 x(n +2 written as in (3.21). Applying (3.41) to this 

representation and performing the integrals over A* via the delta functions, we find 


k 71/ I 2 k 

n / d-'^PK exp(-i ^ plC'Ljgaj) = 5^^(C" • ^). (3.43) 

K=l'^ j=l L=1 

^®The quotient Y always corresponds to the product P 1 P 2 from the gluing procedure (3.18) and thus the 
factor Y{1 — Y)~^ always cancels a factor of [/3i/32(l —/3i/l2)]~^ that arise when the consecutive minors are 
translated into edge variables Oi, /3i and P 2 and that is not present in (3.18). 
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Hence, we can write (3.38) as 


_ f 

d^n+l d^n+2 J 


^^kx{n+2)(jl 

Vol[GL(A:)] 


^n,k ■ W) , 


(3.44) 


where Vl^^k is given in (3.39). 

It would be interesting to further investigate the structure of this expression; we leave 
this for future work. Instead, we will now transform the Grafimannian integral to mo¬ 
mentum twistor space, which will in particular facilitate the explicit calculation of some 
example form factors. 


Momentum twistor space 

Next, we transform our result to momentum twistor space following the strategy of [34, 59]. 



Figure 14: Momenta and dual momenta in the case of form factors, shown for n = 4. In 
contrast to the case of amplitudes, the n on-shell momenta do not add up to zero but to 
the off-shell momentum q of the operator. Hence, the contour is not closed but periodic. 
In order to obtain a closed contour, two on-shell momenta Pn+i and Pn +2 can be inserted 
between any yi and Pi+i of the periodic contour. Two different choices are shown in shaded 
frames. 

In order to introduce the momentum twistor variables Zi = (Aj,^j, 2 i) [60] correspond¬ 
ing to our choice of variables (3.10) and (3.11), we define the dual (super) momenta pi (??i) 
via20 

AiA* = yi- Vi+l ; hUi = '&i- '^i +1 ■ (3-45) 

Note that we base the dual (super) momenta on the closed contour obtained by adding 

Pn+i and Pn +2 instead of the periodic contour as done in [7, 18]; cf. figure 14. Then, we 

define and 2* via the incidence relations 

m = XiVi = km+i, Ui = h'&i = h'&i+i ■ (3.46) 

use the definitions of [59], which coincide with the ones of [7] but differ from the ones of [18] by a 

global sign and a cyclic relabelling. 
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Inverting these relations, we have 


{i+li)gi-i + {ii-l)gi+i + {i-li+l)gi 
{i — li){ii+l) 

{i+li)-gi-i + (ii-l)a^+i + (i-li+l)ai 
— i+1) 


(3.47) 


We start the transformation of the Grafimannian integral from momentum space using 
the representation of as (3.21). We can use part of the GL{k) redundancy 

to fix 


0 ••• 010 
0 ••• 0 0 1 


(3.48) 


As a consequence, the delta functions in (3.21) fix the last two rows of C as 


n! _ \ 1 
L'fc-l i — Ai ) 


G' 

'-'ki 



(3.49) 


Then, (3.38) becomes 




^{k-2)x(n+2) Q/ 

Yol[GL{k - 2) X Tk- 2 ] 


^n,k 52 x(fc- 2 )(^/ . ^4x(fc-2)(^/ . ^ 

(3.50) 


where the integral and the delta functions contain only the first k —2 rows of G'. The shift 
symmetry Tj ._2 acts on these k — 2 rows as 


C'li G'n + n/Ai + ^21^ , I — l,...,/c — 2, 


(3.51) 


with ri/, r 2 i arbitrary. According to [59], (3.47) leads to 


n+2 

E c'uh 


2 = 1 


n+2 

2=1 


n+2 

E c'um 

2=1 


n+2 

/= l,...,A:-2, (3.52) 

2 = 1 


where the matrix D is defined via 


Dii 




(3.53) 


Next, we rewrite the minors of C in terms of minors of D. In [59], it was found that the 
consecutive minors are related as 


{G[...G'^,) = -{l2)---{k-lk){D2... Du-i) (3.54) 

and its natural extension via cyclic shifts. However, we do also need non-consecutive minors, 
as can be seen from (3.39). For these, we find 

{G[... G+++i) = -(12) • • • {k-2 k-l){k-lk+l){D2 ... Du-i) 

-(12) • • • {k-2k-l){k k+l){D2 ... Dk-2Dk) , 

(C[C^ .. C+i) =-(1 3)(34) • • • (A: A:+1)(D3 ... A) 

-{l2){2,A)---{kk+l){D2Di...Dk). 
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Using (3.54), the product of consecutive minors in (3.39) becomes 


( 1 ... k)c> • • • (n+2 • • • k-l)a = 2) • • • (n+2 1))"-'(1 • • • A:)^ • • • (n+2 

For Y, we find using (3.55) 

(n —A:+2 ■ • • n n+l)c'{n+2 1 ■ ■ ■ k — l)c' 


k—l)D ■ 

(3.56) 


Y = 


{n — k+2---n n+2)c'{n+l 1 ■ ■ ■ k — \)ci 

{nn+l){n — k+2> ■ ■ ■n)D 

{nn+2){n—k+2> ■ ■ ■ n)£) + {n+ln+2){n—k+2> ■ ■ ■ n — ln+l)j:) 
{n+2l){l---k-2)D 


(3.57) 


(n+11)(1 • • • k—2)£) + (n+1 n+2) (n+2 2 ■ ■ ■ k — 2)£) 

The remaining steps in the derivation of [59] go through unchanged. First, we use the 
Tk -2 shift symmetry to set Cj^ = Cj 2 = 0, which changes the measure as 


^{k-2)x{n+2) (jl 


= ( 12 ) 


k-2 


^{k-2)x{n) 


(3.58) 


Vol[GL(A; - 2) X Tk- 2 ] Yol[GL{k - 2)] 

Then, we perform the change of integration variables from C' to D, which yields 
Sk-2)x(n)(jf _ /(12)... (n + 21)\^"^ d(fc-2)x(n)^ 

Vol[GL(A: - 2)1 “ V (12p / Vol[GL(A: - 2)] ' ^ 

Finally, we undo the gauge fixing of the first two columns of the C matrix, which yields 
factors of 

{12)S\Dji^) (3.60) 

for I = 1,..., /c — 2. See [59] for details of these steps. 

The final expression we find is 

f cl(*:-2)x(n+2)^ 


Tr 


ri.,2 


Vol[GL(A: - 2)] 


^n,k 


(3.61) 


where 




(n l)(n+l n+2) 


y(i-y) 


-1 


(n n+1)(n+2 1) (1 • • • k-2){2 ■ ■ ■ k-1) ■ ■ ■ {n ■ ■ ■ k-5){n+l ■ ■ ■ k-4:){n+2 ■ ■ ■ k-3) 

(3.62) 

and Y is given in (3.57). 

A convenient choice of reference spinors 

It turns out that one choice of reference spinors is particularly convenient. If we set 

= An+l — -^1 ) Cb — An+2 — ; (3.63) 

the above momentum twistor Grafimannian integral for N^MHV becomes 
^(fc- 2 )x(n+ 2 )^ _y(l_y)-l §Hk-2)\A{k-2) _ 2) 


-F, 


n,2 


Vol[GL(A: — 2)] (1 • • • k — 2) ■ ■ ■ {n - ■ ■ /c —5)(n+l • • • A:—4)(n+2 • • • A: —3) 


with 


Y = 


(n —A;+3 • • • n)(l • • • fc—2) 

(n —A:+3---n—1 n+1)(n+2 2 - ■ ■k — 2) 


(3.64) 


(3.65) 
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Examples at MHV, NMHV and NNMHV 

Let us look at some special cases of the above GraBmannian integral representation. For 
k = 2, the matrix D is zero-dimensional and all consecutive minors of D are 1 whereas all 
non-consecutive minors are 0. Hence, the integral in (3.64) is zero-dimensional while the 
integrand is 1.^^ Considering the prefactor Tn, 2 , this is precisely the correct result. 

For fc = 3, 


D = 


(^di d2 ■ ■ ■ dn+2^ ■ 


(3.66) 


The consecutive minors of D are equal to the single di included in them and the non- 
consecutive minors are equal to the di that is alone on its side of the gap, cf. (3.55). Hence, 
the GraBmannian integral (3.64) becomes 




1 


1 


1 


n,2 


jlx(n-|-2) jj 

Vol[GL(l)] 1 _ d.n+ldr,+2 di---dn dn+ldn+2 




(3.67) 


Note that in all of the examples considered in this subsection we will use the convenient 
choice of reference spinors (3.63) to obtain compact expressions. We have explicitly checked 
that our results are indeed independent of this choice. 

The simplest example for k = 3 is n = 3: 

dlx5^ 


-^3,: 


1 


1 


Vol[GL(l)] 1 - did 2 d 3 d 4 ^d 5 


5^^‘^{di2i -|- d222 -|- d^Z^ -|- ^4.24 -|- d^Z^) . (3.68) 


We can use the GL{1) redundancy to fix ds = (12 34), where the four-bracket is defined as 

{ijkl) = det{ZiZjZi^Zi) = eABCoZf'Z^ Z^ Zj^ (3.69) 

with Zi = (Aj,^j). The remaining four integration variables are then completely determined 
by the delta function: 


Thus, 


di = (2345), d 2 = (3451), d 3 = (4512), d 4 = (5123). 

[12345] 


-^3,3 = -7^3,2 


, (5123)(1234) ’ 

(2 3 4 5)(4 512) 


where the five-bracket is defined as 


[ijklm] = yrx 


5^{{ijkl)r]m + cyclic) 


{ijkl){jkl m){k I mi){I mi j){m i j k) 


(3.70) 

(3.71) 

(3.72) 


This result numerically agrees with the one found in [18]. 
For general n, the denominator of (3.67) has poles for 


di = 0, i = 2 ,...,n —l,n-|-l,n-|- 2 , 


di = 


dn+ldn+2 

dn 


dr), — 




dn+l — 


djdn 

dn+2 


dn+2 — 


didn 

dn+l 


(3.73) 


Note that, although Y is singular when inserting the above values for consecutive and non-consecutive 

— (n—fc-|-3-•-n.)(1-• • fc—2) • -t • .-i • 

———— - - - IS 1 in this case. 


minors, the ratio 


(n—fc-|-3-•-n—1 rH-l)('n+2 2) —(n—fc-|-3---n)(l---fc—2) 
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In principle, one can consider (composite) residues of (3.67) for zero and non-zero values 
of the di- However, we find that it is sufficient to consider residues which are composed of 
individual residues taken at zero.^^ As in the amplitude case discussed in [59], these can 
be characterised by the five dj’s with respect to which no residues are taken. In contrast to 
the amplitude case, these have to include di and dn- We have to consider two cases. In the 
hrst case, no residues are taken with respect to dn+i and dn+ 2 - The resulting expressions 
are 


Resj = 


1 


1 (rH-2 lni){lni 

{n i rtfl rjH-2) {i rH-1 rt\^ 1) 


[zn+ln+21n] 


(3.74) 


where i G {2,... ,n — 1}. In the second case, at least one residue is taken with respect to 
either dn+i or dn+ 2 - The resulting expressions are 


R-esij,fc = [ijkln] 


(3.75) 


where i, j, /c G { 2 ,... , n — 1 , n + 1 , n + 2 }. 

An additional property of form factor top-cell diagrams arising here is that we have to 
take the sum of more than one form. This can also be achieved by shifting the legs between 
which the minimal form factor is glued-in from (re, 1 ) to (re -|- s mod re, 1 -|- s mod re), cf. 
hgure 14. 

Numerically comparing with the results of [18], we find 
J- 4^3 = T”4,2(+Res3 -|- Res2,3,5 -|- Res 3 ^ -|- Res2_3^5), 

— T5j2(+Rgs 4 -|- Res3^4^6 + ^ 682 ^ 3^6 T RGS 3 ^ — Res2,3,4 (3.76) 

-|- Res2,3,6 + RGS 3 4 7 — Res2^3^4 + Res 5 ^) , 

where the superscript s specifies the shift. This also gives further support for the relation 
(3.3). 

Finally, we look at the simplest example of A: = 4, namely re = 4. In this case, the 
matrix D can be gauge-fixed to be 


D = 


/l 0 di3 dii di5 die 

yO 1 ^23 ^25 d26 


The delta functions completely fix their entries to 

^ (i456) ^ (i356) (i346) , (i345) 

“ “(3456) ’ “ “^(3456) ’ “ “(3456) ’ “ “^(3456) 

where i = 1,2. Hence, 

(1345)(1346)(1356)(2346)(2356)(2456) [13456] [23456] 
(1234)(1236)(3456)2((1246)(1345) + (1 256)(3456)) 


(3.77) 


(3.78) 


(3.79) 


We have successfully checked (components of this expression) against [7]. 

^^This can also be understood from the corresponding on-shell diagrams. 

^^Note that there is a subtle sign occurring in the evaluation of the momentum twistor Grafimannian 
integral for J- 4 +, which is also present in the case of the corresponding amplitude He, 4 . 
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4 Integrability and form factors 

After having already used the integrability-related R operators in section 2, we now study 
the integrable structure of form factors in A/" = 4 SYM theory more carefully. Apart from 
linking integrability approaches from the spectral problem and the study of amplitudes, 
this is also motivated by the search for symmetries. 

We approach this problem by introducing the spin-chain monodromy matrix as it 
appeared in the context of tree-level amplitudes [35, 36]. While the on-shell part of the form 
factors, studied in the previous sections and built from the R operators (2.14), is Yangian 
invariant, we find that this symmetry is broken by the insertion of the minimal form factor, 
i.e. the off-shell part. However, the off-shell part can be interpreted as an eigenvector of the 
corresponding transfer matrix. From this we show that all form factors of the chiral stress- 
tensor super multiplet are annihilated by the corresponding transfer matrix.This is the 
analogue of Yangian invariance for form factors. Finally, we show that the transfer matrix 
acts diagonally on a given minimal form factor of a generic operator if the corresponding 
operator renormalises multiplicatively, i.e. is an eigenstate of the integrable Hamiltonian 
studied at one-loop order in the spectral problem. As a consequence, also certain planar 
leading singularities of loop-level form factors of generic operators are eigenstates of the 
transfer matrix. 


4.1 Spin chains and Yangian invariance 


In the following, we introduce the spin chain that appeared in the context of tree-level 
amplitudes [35, 36, 39, 40]. In the spin-chain language, the integrability construction is 
naturally formulated using the complex Lie algebra g[(4|4) instead of psu(2,2|4). 

The 0[(4|4)-invariant Lax operator relevant for the construction of the integrable spin 
chain naturally acts on the tensor product of two spaces and depends on the spectral 
parameter u: 

£i{u)=--\-- + xfpf, (4.1) 

i 


where A = (a, d, A). We have introduced a graphical notation, usually used in the context 
of vertex models, to depict the tensor structure of the Lax operator, see e.g. [36], and | • | 
denotes the grading. While the auxiliary space is finite-dimensional with the generators 
(e'^^)cD = and illustrated by the dashed horizontal line, the quantum space at site 

i is infinite-dimensional and is denoted by the vertical line. The corresponding generators 
are realised in the Jordan-Schwinger form using the Heisenberg pairs 


X' 


A 


d 9 \ 

’ (9A“ ’ dfj^ ) ’ 


p. 


A 



with [x-^,p^] = (-l)l-^l(I-^®, 

(4.2) 


where [•, •] denotes the graded commutator; see e.g. [37]. 

The spin-chain monodromy matrix is built from the n-fold tensor product of the Lax 
operators (4.1) in the infinite-dimensional quantum space and matrix multiplication in the 


^■^For this particular super multiplet, the corresponding eigenvalue is zero. 
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auxiliary space. Graphically, multiplication from the left or the right in the auxiliary space 
(quantum space) correspond to attaching vertices from left (bottom) or right (top). We 
define 

Mn{u, {Ui}) = = £„(u - Vn) ' ■ ■ C2{u - V2)Ci{u - Ui) , 

n 2 1 

with inhomogeneities Vi that are local shifts of the spectral parameter u. 

For later purposes, we also introduce the corresponding transfer matrix, which is con¬ 
structed from the monodromy matrix as the super trace over the auxiliary space: 

Tn{u, {uj}) = = str7W„(u, {ui}). (4.4) 

n 2 1 

As a consequence of the Yang-Baxter equation, this transfer matrix is g[(4|4) invariant: 

n 

[nu,{v,}),J2j^^]=0. (4.5) 

i=l 

Tree-level scattering amplitudes are Yangian invariant [43]. Instead of using Drinfeld’s 
first realisation, this can be expressed as a set of eigenvalue equations involving the mon¬ 
odromy matrix in (4.3), cf. [35, 36]: 

M.n{u,{vi))A(x.lLA. (4.6) 

While for the physical amplitude the inhomogeneities Vi are zero, they can be set to non¬ 
zero values to obtain deformations of the amplitude with non-vanishing central charges 
[37, 55]. 

From (4.6), it follows that the transfer matrix (4.4) acts diagonally on tree-level am¬ 
plitudes yielding a vanishing eigenvalue. However, this condition is weaker than the set of 
eigenvalue equations in (4.6). 

4.2 Form factors of the chiral stress-tensor multiplet 

Next, we study the action of the monodromy matrix (4.3) on more general form factor 
expressions of the chiral stress-tensor multiplet, which we define as 

with 

:Fl,{k J^2,2{k - 1, k) 5^^, • • • 5- . (4.8) 

They are constructed from a chain of R operators acting on a vacuum state that is composed 
of the amplitude vacua 5~ as well as the minimal form factor (2.20). The R operators 
and amplitude vacua correspond to the on-shell part of the diagram, with the minimal form 
factor cut out. The number m of R operators depends on the diagram under consideration. 
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Moreover, the R operators have to be chosen in such a way that the corresponding diagram 
is planar.This generalises the construction we used in section 2 (see e.g. (2.18) and 
figures 4, 5, 9), and is another way of writing the gluing of diagrams that was performed in 
section 3. Note that these objects correspond to any planar on-shell diagram containing an 
insertion of the minimal form factor, including top-cell diagrams, individual BCFW terms, 
factorisation channels etc. 


On-shell part 

It was discussed in [35, 39, 40] that tree-level scattering amplitudes can be constructed via 
the method of R operators which naturally include the inhomogeneities Vi as deformations 
of the local central charges [37].^® These operators, defined by their action on functions of 
the kinematic data in (2.14), can be formally written as 

Ry(u)= (49) 

j i 

The operator R can be seen as one of two basic building blocks for Yangian invariants. It 
satishes the Yang-Baxter equation 


Ui) , (4.10) 

with the Lax operators defined in (4.1). This equation can be depicted as 

(4.11) 

j i 3 i 



The other basic building blocks are the vacuum states 5f = (f^(Ai) and 5^ = (f^(Ai)(5^(r)j) 
introduced in (2.15). They satisfy 


Ci{u) 5f = (u - 1) 1 5f , Ci{u) = n 1 , 


which can be depicted as 



+ 

- 1) I , 

i 



(4.12) 


(4.13) 


The properties (4.10) and (4.12) guarantee that an appropriate combination of R operators 
with a suitable choice of inhomogeneities acting on the tensor product of vacuum states 
6f and 5'^ is Yangian invariant [35, 39, 40]; the required choice of inhomogeneities will be 

^®In particular, we assume that for each operator Rij the indices satisfy i < j', this corresponds to the 
chosen convention for the BCFW shifts. 

^®See also [36], where Yangian invariants were constructed using Bethe ansatz methods. 
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discussed in the following. However, further below we will also see that Yangian invariance 
is broken by the insertion of the minimal form factor. 

As discussed in [35, 39, 40], the monodromy matrix (4.3) satishes certain exchange 
relations with a chain of R operators. As a consequence of (4.10), one finds 

M(u, {ui}) Rqji(^i) ■ • • Rimim(^m) ~ (^i)''' ^imjmi^m) A4(tt, , (4.14) 

where A4(u, {^^-(i)}) denotes the monodromy matrix in (4.3) with the inhomogeneity at 
site i permuted such that Vi is replaced by Uo-(j). Here, a is the permutation associated to 
the on-shell diagram and can be read off as discussed in section 2.1. The inhomogeneities 
Vi and the spectral parameters Zi have to satisfy the relations 


Zi = Vrtiii) - with Tf, = (H, ji) <]•••<] (if, j^), £ = 1,..., m , (4.15) 

see [35, 39, 40].^^ The inhomogeneities Vi associated to the i*'^ external leg, i.e. to site i, 
are related to the central charges c* via [38] 

a = Vi- . (4.16) 

Therefore, for a planar on-shell diagram with valid deformations, we can commute the 
monodromy matrix through the chain of R operators in (4.7) using (4.14): 


Mniu,{Vi})J^ = Riijiizi) ■ --Ri^j^iZm) A4„(u, }) J'2,2(^ “ 1,^) • (4-17) 

Since the Lax operators act diagonally on the vacua (4.12), we can eliminate all Lax 
operators that do not act on the minimal form factor. We end up with the monodromy 
matrix 

M2{u, {u^(j)}) = Ck{u - u^(fc))£fc_i(u - u^(fc_i)) (4.18) 

of length two acting on the minimal form factor J- 2 , 2 {k — 1, k) at sites k — 1 and k: 
Mn{u,{v„^i)}) J^l2ik - l,k) 




' ^k-2 


(4.19) 


^k+l 


^(i)) • 


(4.20) 


M2{u, {v„(i)})J^2,2{k - 1, k) 

where 

k—2 n 

f(^^,{Mo}) = “1) n 

i=l j=k+l 

As we will discuss below, this significant difference to tree-level amplitudes breaks the 
Yangian invariance. However, as we will also see, some of the integrable structure remains. 

Using the graphical language introduced earlier, we depict the formulas discussed above 
for the case of J-n ,2 in figure 15. The left picture in figure 15 represents the monodromy 
matrix (4.3) acting on the chain of R operators (BCFW bridges) as introduced in (2.14) 
and (4.9) contracted with the minimal form factor (2.3) and the corresponding vacua (2.15). 
In (4.17), we commuted the monodromy matrix through the chain of R operators (BCFW 
bridges) as shown in the middle picture. The action of the Lax operators on the delta 
functions of the on-shell vacua was discussed in (4.12). As in (4.19), we end up with a 
monodromy matrix of length two acting only on the minimal form factor as shown in the 
right picture of figure 15. 


^^Recall that m denotes the number of R operators. 
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Figure 15: Action of the monodromy on 


Minimal form factor 

One can explicitly check that the minimal from factor is not an eigenstate of the length- 
two monodromy matrix (4.3), and thus not Yangian invariant, cf. (4.6). This can be seen 
since, for example, the momentum-like generators do not contain the off-shell momentum 
q. Acting on the minimal form factor, this produces 

+ ^k^k)^^{^k-l^k-l + ^k^k ~ <?) j (4-21) 

which does not vanish, as would be required for Yangian invariance in the sense of (4.6). 
However, as we will show, the minimal form factor is annihilated by the graded sum of 
the Yangian generators on the diagonal of the monodromy matrix, i.e. the transfer matrix 
(4.4), for equal inhomogeneities: 

T2{u - v) = StT Ckiu - V„(^k))Ck-l{u - , with = Va{k) = V ■ (4.22) 

This can be seen as follows. First, note that we can look at the action of the transfer 
matrix (4.22) on a single component of J^2,2 due to its g[(4|4) invariance (4.5). We take as 
an example the component 




(4.23) 


which corresponds to the minimal form factor of the operator in (1.2) with 

outgoing scalars label hy k — 1 and k. Note that the transfer matrix does not depend on 
7“. Second, one can check that the transfer matrix annihilates the momentum-conserving 
delta functions, 

T2{u-v)6^{XiXi + X2X2-q) =0, (4.24) 

and thus acts only on the product of fj's, yielding 


T2{u - 

= ((u - u - l)(x^_ip^_i + xfp^) + 


?;+l _ Q 


Vk-iVk-iVk % 


(4.25) 
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This shows that the minimal form factor is an eigenstate of the transfer matrix with 
eigenvalue zero, i.e. 

r2(u-u)J2,2 = 0. (4.26) 

Moreover, due to the fact that the monodromy matrix, and therefore also the transfer 
matrix, commutes through the chain of R operators as discussed above (4.14), the same 
statement applies to any planar on-shell diagram with an insertion of the minimal form 
factor (4.7): 

Tn{u,{vi})P = (4.27) 

with the constraints on the inhomogeneities given in (4.22). As for amplitudes, the whole 
argument is in particular valid for vanishing inhomogeneities, i.e. for undeformed form 
factors: 

TnJ^n,k = 0 . (4.28) 

Note the similarity to the Yangian invariance condition of scattering amplitudes (4.6): 
taking the super trace for them shows that they are also eigenstates of the transfer matrix 
with eigenvalue zero. Equation (4.27) implies that form factors, although not Yangian 
invariant, are still annihilated by a certain combination of the Yangian generators. 


4.3 Generic operators 


In the previous section, we have shown that form factors of the stress-tensor super multiplet 
can be interpreted as eigenvectors of the transfer matrix (4.4) with vanishing eigenvalue. 
In the following, we want to extend this to all operators. We will study the action of the 
homogeneous transfer matrix on the minimal form factors of generic single-trace operators, 
and find that they too are eigenvectors provided that the single-trace operators are chosen 
as eigenvectors of the spectral problem. Combining this with an R-operator construction 
similar to the one in the previous section, it follows that all planar on-shell diagrams that 
include these minimal form factors are also eigenstates of the transfer matrix. In general, 
these objects should correspond to leading singularities of loop-level form factors. For the 
purposes of the following discussion it will be convenient to work in components and, in 
contrast to the rest of this article, we do not employ harmonic superspace variables here. 

Generic operators in AA = 4 SYM theory can be conveniently represented via two sets of 
bosonic oscillators and and one set of fermionic oscillators acting on a suitable 
vacuum [61, 62]. The oscillators in this oscillator picture transform under psu(2,2|4) in the 
same way as the super spinor helicity variables Xf , Xf and fjf and the algebras are formally 
the same if one identifies 






bf Af, 




r. _ d P, _ d 


dX? 


dXf ’ 


M,A 


di,A — 


d 

dvf ’ 


(4.29) 


see [63] for a detailed comparison of these representations. 

In [64], this identification was used to connect the one-loop dilatation operator to the 
tree-level four-point scattering amplitude based on symmetry considerations. The field- 
theoretic quantities behind such an identification in the composite operators are actually 


36 





form factors. Concretely, it was shown in [19] via an explicit Feynman diagram calcula¬ 
tion that the colour-ordered minimal tree-level super form factors of generic single-trace 
operators O can be obtained from their representation in the oscillator picture as 


■■■ ,L]q) 
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o 
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(4.30) 


where L is the number of fields in the single-trace operator which has been translated 
according to (4.29). 

Due to the g[(4|4)-invariance of the transfer matrix (4.5), it commutes with any 
function /(X)ii=i This in particular implies that it commutes with the momentum- 

conserving delta function in (4.30), which becomes clear after rewriting 








(4.31) 


and recalling that = Af Af. Thus, the transfer matrix only acts on the operator trans¬ 
lated into spinor helicity variables. Translating the spinor helicity variables in the transfer 
matrix T to oscillators using (4.29) yields 


T(tt) = str L2,(ri) ■ • • Li(u) with Li(u) = £j(n) 


A? 








9i,A, vf 


We obtain the relation 


where 


Tl{u)Fo,L = J^Tl{u)0,L 
-T Tl{u)o,l = L5‘^ AjAj — q 
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(4.32) 

(4.33) 

(4.34) 


Having expressed the action of the transfer matrix T on the minimal form factor J^o,L 
in terms of the transfer matrix Tl, cf- (4.33), we will now argue that acts diagonally 
on O, i.e. 

Tl{u) 0 = t{u)0, (4.35) 


if the state corresponding to the operator O is an eigenstate of the spin-chain Hamiltonian 
H, i.e. the one-loop dilatation operator of A/" = 4 SYM theory. Here, t{u) is a polynomial 
in the spectral parameter While the particular transfer matrix T(m) does not contain 
the spin-chain Hamiltonian H, it is commonly used to diagonalise the commuting family of 
operators [65]. Here, we show that for Vi = 0 the transfer matrix T(u) commutes with H 
using a criterion by Sutherland [66] and therefore belongs to the same family of commuting 
operators, see also [67]. This is a consequence of the Yang-Baxter equation 


Rj,i+i(tt)Lj(u -h u')'Li+i{u') = Li+i(u')Li(u -h u')Rj,j+i(u) (4.36) 

^®In the framework of the Bethe ansatz, the eigenvalues !(«) can be parametrised by the Bethe roots. 
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studied in [37] where the harmonic R matrix R was derived, see also [36]. The expansion 
of the harmonic R matrix contains the Hamiltonian density T-Li^i+i at first order in the 
spectral parameter = Vi^i+i{l+u'Hi^i+i + ...), see e.g. [67] as well as [37, 68] 

where this relation was discussed in relation to A/" = 4 SYM theory. Taking the derivative 
of (4.36) with respect to u and subsequently multiplying with the permutation operator 
Vi^i+i, we obtain 

[ni,i+i,U{u')U+i{u')] = U{n') - U+i{u') . (4.37) 

As a consequence, one finds that the Hamiltonian commutes with the transfer matrix 

L 

[H, T(u)] = 0, with H = ^ , (4.38) 

i=l 

where periodic boundary conditions 'Hl,l+i = T~iL,i are imposed. 

Just as for the minimal form factor of the chiral stress-tensor multiplet in section 4.2, 
we can glue planar on-shell diagrams to the minimal form factor (4.30) using R operators, cf. 
(4.14). By construction, this part is Yangian invariant and the monodromy matrix can be 
commuted through the chain of R operators as shown in (4.17). The commutation relations 
for the transfer matrix built from that monodromy matrix follow immediately after taking 
the trace in the auxiliary space, see the first step in figure 16. Using (4.33), we find 

Tn{u) Po,n = f{u)^TLO,n , (4.39) 

where f(u) denotes the factors arising from the action of the Lax operators on the vacuum, 
see (4.12). Finally, from the argument presented above, we find that these generalisations 
of (4.7) are eigenstates of the transfer matrix T, 

Tn{u) To,n = = f{'u)i{u) To,n , (4.40) 

if the operator satisfies the eigenvalue equation (4.35). This generalises the corresponding 
identity (4.27) for the stress-tensor super multiplet and can be denoted graphically as 
shown in figure 16.^® 

Note that J-o,n does not necessarily correspond to a tree-level form factor of the com¬ 
posite operator O. However, at least some of the on-shell diagrams correspond to leading 
singularities of loop-level form factors. It would be interesting to see whether general tree- 
level form factors are eigenstates of the transfer matrix, and whether the identities (4.39) 
for leading singularities are hints of similar integrability properties at loop level. 

5 Conclusion and outlook 

In this paper, we have extended many concepts that were developed in the context of the 
purely on-shell amplitudes to the partially off-shell form factors, focussing on the tree-level 
form factors of the chiral part of the stress-tensor multiplet as an example. 

^®In fact, the intermediate state in figure 16 depicts the generalisation of the right hand side of (4.19) to 
the transfer matrix, which coincides with the intermediate step in (4.40) via (4.33). 
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Figure 16: Action of the transfer matrix on an on-shell diagram with an insertion of the 
minimal form factor of the operator O. 


We have shown that on-shell diagrams can be used to characterise form factors by 
including the minimal form factor as a further building block in addition to the three- 
point MHV and MHV amplitudes. Apart from the equivalence moves for amplitudes, this 
requires a rotation move that reflects the cyclic invariance of the three-point form factors. 
Moreover, we can extend the concept of a top-cell diagram to form factors. Whereas one 
top-cell diagram suffices for amplitudes, we require several ones for form factors. We have 
given a conjecture for the top-cell diagrams for all numbers of on-shell particles n and MHV 
degree k, which is based on a relation to the amplitude case. We have explicitly checked 
this conjecture against the known results for all n at MHV level, up to five external points 
in the NMHV sector and for the simplest example at NNMHV. 

Moreover, we have rewritten the previously obtained expressions into the form of a 
(deformed) GraBmannian integral. As we use two on-shell momenta to parametrise the 
off-shell momentum of the composite operator, it in general involves the GraBmannian 
G{k, n+2), where n is the number of external on-shell states. This construction geometrises 
the (super) momentum conservation and is based on gluing the minimal form factor to the 
rest of the on-shell diagram. We have obtained the GraBmannian integral in spinor helicity 
as well as twistor and momentum twistor variables. As can be seen from (3.39) and (3.38), 
there are significant differences between the GraBmannian integral formulas for scattering 
amplitudes and form factors. The form in the case of scattering amplitudes is expressed in 
terms of consecutive cyclic minors labelled by the external particles. The gluing procedure 
for form factors, however, results in a form that contains minors involving non-consecutive 
labels as well. Moreover, we use two fictitious on-shell particles to parametrise the off-shell 
momentum of the operator and hence the minors of the form factor GraBmannian are 
labelled not just by physical external particles but by two additional ones and it seems 
generally not possible to disentangle these ones from the rest. 

Introducing a central-charge deformation to form factors, we could construct them via 
R operators in analogy to the amplitude case. While amplitudes are Yangian invariant and 
hence eigenvectors of the monodromy matrix, the behaviour of general n-point form factors 
when acting with the monodromy matrix is entirely determined by its residual action on 
the minimal form factor and hence the corresponding Yangian, see figure 15. In order to 
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obtain an eigenvalue equation for the form factor, we have taken the super trace of the 
homogeneous monodromy matrix, which yields the homogeneous transfer matrix. In partic¬ 
ular, we find that the minimal form factor of the stress-tensor super multiplet is annihilated 
by this transfer matrix, which contains a subset of the Yangian generators. This equation 
is shown in (4.28). Our construction of the n-point tree-level form factors of the chiral part 
of the stress-tensor multiplet via the integrability-inspired method of R operators shows 
that they satisfy (4.28) as well. Furthermore, we have shown that the minimal form factors 
of all operators can be interpreted as eigenstates of the homogeneous transfer matrix and 
established a connection between the integrable structure of the spin chain that appeared 
in the spectral problem and the one that appeared in the study of tree-level amplitudes. 
Finally, we have argued that the minimal form factor of generic operators can always be 
dressed with a chain of R operators without spoiling the eigenvalue equation with respect 
to the transfer matrix, see (4.39). The resulting objects can be interpreted as leading sin¬ 
gularities of loop-level form factors. 

The results discussed in this work open up very interesting directions of future research. 
Our construction of the form factor top-cell diagram uses the corresponding diagram from 
scattering amplitudes with a box replaced with the minimal form factor, see (3.3). This 
construction has worked in all our examples and can be proven for individual BCFW terms, 
but it is desirable to have a general proof that the conjectures (3.3) and (3.38), (3.39) for 
the top-cell diagrams and the corresponding Grafimannian integrals for all N^MHV form 
factors indeed produce all possible BCFW terms. 

As mentioned earlier, we have tested our Grafimannian formula with known results 
based on case studies. For the examples we have checked, we could determine the combi¬ 
nation of the residues which gives the correct tree-level form factors, but it would be very 
interesting to find a general prediction of a contour which gives the right combination. For 
the scattering amplitudes, such a contour was determined by the twistor string theory for¬ 
mulation of scattering amplitudes [69-71], which leads to the question whether form factors 
also have an interpretation in terms of an underlying twistor string theory. The presence of 
non-consecutive minors as well as the fact that we have multiple top-cell diagrams makes 
the classification of residues a more formidable problem. 

So far, a deeper understanding of the geometry of the Grafimannian formulation is 
missing. In particular, it would be interesting to see whether the form in (3.38) follows 
from a modified notion of positivity. 

Another fruitful direction to pursue would be to extend our result to loop level and 
to obtain a “formfactorhedron” as a counterpart to the amplituhedron [72]. It is known 
that even for planar form factors at two-loop level we need Feynman diagrams that are 
non-planar after removing the minimal form factor. It will be interesting to see how such 
apparent non-planarity for loop-level form factor plays a role in the on-shell diagrams and 
Grafimannian formulation.^*^ Moreover, there has been very interesting progress in study- 

®°For a discussion of non-planar on-shell diagrams for amplitudes, see [73-77]. In particular, non- 
consecutive minors also appear in this context. 


40 



ing loop-level correlation functions in A/” = 4 SYM theory using the Lagrangian insertion 
techniques [78] and it may be interesting to understand a similar picture for form factors at 
loop level, or to investigate whether the minimal form factor insertions even have a direct 
interpretation within this framework. 

While we have focused on the form factor of the chiral part of the stress-tensor super 
multiplet in first half of this paper, we are convinced that our results can be extended 
to general operators; in particular to operators that are non-protected and have a length 
L > 3. 

We expect that the minimal form factor allows for deformations consistent with the 
deformation introduced through the R operators. The resulting eigenvalue equation for the 
inhomogeneous transfer matrix suggest that any given on-shell diagram with an insertion 
of a deformed minimal form factor is fully characterised by the inhomogeneities and Bethe 
roots and as such can be determined using Bethe ansatz methods along the lines of [36]. 
Hopefully, this will yield a uniform description of the observables of planar A/" = 4 SYM 
theory as an integrable model at weak coupling and beyond, see also the recent and very 
promising approach in [79, 80]. Constructing the form factors as solutions to (4.39) at 
loop level should in particular yield an integrability-based approach to the eigenvalues 
and eigenvectors of the dilatation operator.While the eigenvalues can be obtained to 
very high loop orders via integrability, the corresponding eigenvectors are only known to 
one-loop order for generic operators. 

Going beyond form factors, it would be very interesting to glue the on-shell diagrams 
of form factors together to obtain on-shell diagrams for correlations functions, which are 
purely off-shell objects. The results of [13] suggest that this is meaningful at least at the 
level of leading singularities. 

Finally, it would be interesting to extend our results to other theories. In particular, 
both form factors [82-85] as well as (deformed) on-shell diagrams and Grafimannians [41, 
86, 87] were already studied in ABJM theory [88] and could be combined as was done in 
this paper for J\f = A SYM theory. 
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